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NOTATION 


significant  wave  amplitude 


Charnock  proportionality  constant 
Wave  phase  velocity 

T 


wind-stress  coefficient 


Froude  number 


Fetch 

Gravitational  acceleration 
Significant  wave  height 


Dynamic  roughness 
=  30  r\,  surface  roughness 

,  shear  velocity 

Wind  velocity 

Wind  velocity  measured  at  distance  y  above  mean 
water  level 

=  0.42,  Karman  universal  constant 
Wave  length 
Density  of  air 

Surface  tension  at  air-water  interface 
Wind  stress 
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ABSTRACT 


Based  on  the  compiled  data  of  thirty  independent  oceanic 
observations,  this  report  systematically  presents  the  wind- 
stress  coefficient,  the  surface  roughness  and  the  boundary 
layer  flow  regime  at  the  air-sea  interface  under  various  wind 
conditions.  The  air  flow  near  the  water  surface  is  shown  to  be 
aerodynamically  rough  or  in  transition  region  except  at  very 
low  wind  velocity  (Uio  <  3  m/sec).  Both  the  wind -stress  coef¬ 
ficient  and  the  surface  roughness  are  found  to  increase  with 
the  wind  velocity  when  Uio  is  less  than  IB  m/sec  and  reach  a 
saturated  value  for  Uio  greater  than  15  m/sec.  Based  on  the 
oceanic  wave  observations,  it  shows  that  the  presence  of  this 
discontinuity  at  Uio  =  15  m/sec  is  due  to  an  increase  in  the 
wind  velocity  (measured  at  the  significant  wave  amplitude  above 
the  mean  water  level)  beyond  the  average  wave  phase  velocity. 
This  finding  provides  a  well-defined  separation  for  the  often 
quoted  terms  "light"  and  "strong"  winds,  and  also  explains  the 
existence  of  the  so-called  "critical  wind  velocity"  suggested 
by  Munk.  It  also  shows  that  the  surface  roughness  is  governed 
by  the  amplitude  of  the  short  gravity  waves  rather  than  the 
mean  square  surface  slope.  Charnock's  relationship  is  shown  to 
be  applicable  to  most  of  the  oceanic  data, and  Charnock's  pro¬ 
portionality  constant  is  determined, 

— - —  =  0.0156. 


u*/  g 
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Finally,  two  approximate  formulae  for  the  wind-stress  coeffi¬ 
cient 

i 

Cio  =  0.5  x  uf0  x  10" 3  for  light  wind 


and 


Cio  =  2.6  x  10“ 3  for  strong  wind 


are  suggested  for  oceanic  applications. 


INTRODUCTION 


The  wind-stress  acting  at  the  air-sea  interface  locally  is 
an  important  influence  on  the  wind-wave  interaction,  including 
the  generation  of  water  surface  set-up,  drift  current  and  sur¬ 
face  waves,  as  well  as  on  the  transfer  of  heat  and  mass  through 
the  Interface.  For  the  global  scale,  therefore,  it  influences 
the  geophysical  state  of  the  atmosphere.  Numerous  experimental 
investigations,  analytical  studies,  and  summarized  reviews  of 
this  subject  have  appeared  in  various  journals,  showing  the 
great  effort  that  has  been  made  in  this  area,  and  indicating 
the  Importance  of  the  subject. 

In  considering  the  difficulties  involved  in  field  wind- 
stress  determinations,  discrepancies  among  their  results  are 
Inevitable.  However,  the  trend  of  the  data,  hidden  by  the  dis¬ 
crepancies,  may  be  uncovered  if  many  sets  of  data  are  compiled, 
and  if  the  compiled  data  are  analyzed  with  the  help  of  physical 
reasoning.  Extending  a  previous  laboratory  study  (Wu,  196?) 
about  the  physics  of  wind-stress  and  surface  roughness  at  the 
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air-water  Interface,  and  a  subsequent  work  about  the  scaling  of 
the  wind-stress  coefficient  (Wu,  1968),  the  abundant  field  data 
are  herein  correlated  in  order  to  investigate  the  wind-wave 
interactions  at  the  air-sea  interface  and  to  establish  a  reli¬ 
able  formula  or  formulae  for  oceanic  applications. 

The  wind  velocity  near  the  air-sea  interface  was  reported 
unanimously  to  follow  the  logarithmic  distribtion,  see  Roll 
(1965).  This  distribution  has  also  been  verified  by  Wu  (1967) 

In  a  laboratory  channel.  In  the  same  study,  an  expression  re¬ 
lating  surface  roughness  and  wind-stress,  first  suggested  by 
Charnock  (1955)>  has  been  established,  provided  that  the  sur¬ 
face  roughness  is  governed  by  gravity  waves.  Characterizing 
the  equilibrium  condition  between  the  wind  and  waves,  this  non- 
dimensional  expression  coordinates  satisfactorily  the  labora¬ 
tory  and  field  determined  wind-stress  coefficients,  despite  the 
great  differences  In  7:eir  fetches  (Wu,  1968). 

Based  on  the  collected  oceanic  data  and  by  accepting  the 
logarithmic  nature  of  the  wind  profile,  we  have  studied  the 
applicability  of  Charnock's  expression  and  determined  Charnock':; 
proportionality  constant.  The  present  determined  value  ( 0 .  t  1 96 ) 
is  much  smaller  than  the  commonly  accepted  one  (0.078)  previ¬ 
ously  obtained  by  Hay  ( 1951*)  •  (His  measurement  consists  of  six 
wind  profiles  and  contrioutes  six  points  to  the  preset  ’  data 

point  collection).  Moreover,  this  report  systematic'1  1  !.v  pra  - 
sents,  for  the  first  time,  the  wind-stress  coef  n  c ! »  n-  •  \d  the 
surface  roughness  at  the  air-sea  interface  under  difT1  w'.ud 
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condltlons.  Furthermore,  we  have  shown  and  explained  the  exis¬ 
tence  of  a  discontinuity  of  the  wind-stress  coefficient  at 
Uio  =  15  m/sec,  where  U10  is  the  wind  velocity  measured  at  10  m 
above  the  mean  water  level.  We  relate  this  discontinuity  to  the 
"critical  wind  velocity",  which  was  first  suggested  by  Munk  (1947), 
and  has  been  misinterpreted  as  being  due  to  the  transition  of  the 
boundary  layer  from  laminar  to  turbulent.  We  suggest  this  wind 
velocity  (Uio  =  15  m/sec)  as  the  separation  for  the  often  quoted 
but  very  confusing  geophysical  terms  "light"  and  "strong"  winds. 
Based  on  all  of  these  findings  and  physical  reasonings,  two  ap¬ 
proximate  formulae  are  suggested  and  they  are  shown  to  have  much 
less  deviation  from  the  compiled  wind-stress  data  than  the  formu¬ 
lae  proposed  by  other  investigators. 

PREVIOUSLY  PROPOSED  CORRELATION 

Since  the  validity  of  the  logarithmic  wind -veloci ty  distri¬ 
bution  with  height  is  well  established,  the  wind-stress  coeffi¬ 
cient,  C  ,  defined  in  terms  of  the  wind  velocity,  U  ,  measured 

y  y 

at  a  certain  height,  y,  can  be  written  as 


where  p  is  the  air  density,  k  Is  the  Karma n  constant,  and  h  is 
the  virtual  origin  of  the  logarithmic  velocity  profile  (or  the 
dynamic  roughness). 
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It  was  found  in  laboratory  experiments  (Wu,  1967)  that  at 
high  wind  velocities  when  the  wind-stress  is  supported  by  the 
form  drag  of  the  gravity  waves,  a  nondimensiona 1  expression, 
first  suggested  by  Charnock  (1955),  could  be  established 

— - —  =  a  =  0.0112  .  [2] 

u3*  /g 

By  substituting  [2]  into  [ l] ,  an  equation  for  determining  the 
wind-stress  coefficient  was  obtained  (Wu,  1968) 

[31 

in  which  F  is  defined  as 


This  equation  suggests  the  Froude  scaling  law  for  wind-stress 
coefficients  and  has  been  shown  by  Wu  (1967,  1968)  to  correlate 
satisfactorily  tne  results  obtained  from  measurements  at  all 
fetches . 

COMPIT  ATION  OF  WIND-STRESS  DATA 

The  results  of  forty-two  experimental  Investigations,  twelve 
laboratory  studies  and  thirty  oceanic  observations,  have  been 
collected.  Most  of  the  data  were  obtained  from  Francis  ( 1  '51 ,  l-^-'O  * 


— i—  =  -  -in 

vr~  * 

y 
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Hill  (1962),  Phillips  (1966),  Roll  (1965),  and  Wilson  (i960) 
instead  of  the  original  articles;  no  attempt  is  made  here  to  ex¬ 
amine  the  accuracy  of  the  individual  works  and  the  various 
methods  employed  and  equal  weight  has  been  assigned  to  eacn  data 
point . 

The  compiled  data  (C^  vs  U  )  are  plotted  in  Figure  1,  in 
which  different  symbols  are  used  for  the  laboratory  and  tne 
field  results.  For  the  latter,  the  wind  velocity  was  measured 
at  a  unique  height  above  the  mean  water  level,  y  =  10  m.  As  for 
the  former,  six  authors  adopted  y  =  10  cm;  the  rest  of  the  data 
were  presented  with  the  average  wind  velocity  across  the  tunnel 
section.  However,  it  is  believed  that  for  an  ordinary  wind  wave 
channel  these  two  wind  velocities  may  not  differ  to  a  great  ex¬ 
tent  and  therefore  are  used  interchangeably.  This  assumption 
will  certainly  limit  the  accurate  evaluation  of  the  laboratory 
data;  however,  the  main  purpose  of  this  paper  is  to  t-xamine 
oceanic  observations. 

Much  scattering  of  the  laboratory  results  is  shown  in  Fig¬ 
ure  1,  but,  surprisingly  enough,  the  oceanic  data  show  consider¬ 
able  consistency  among  themselves.  These  two  groups  of  data 
(laboratory  measurements  and  oceanic  observations)  are  distinctly 
separated  from  each  otner  largely  due  to  great  difference  in  their 
fetches, and  they  can  be  correlated  by  means  of  [3]  witn  the  sub¬ 
stitution  of  y  =  10  cm  for  laboratory  results  and  y  =  10  m  for 
oceanic  measurements;  two  continuous  lines  are  drawn  in  Figure  1 
to  show  this  correlation.  For  readers  Interested  in  scaling  the 
wind-stress  coefficient,  a  detailed  discussion  is  offered  oy 
Wu  (1968). 


X, 
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The  oceanic  data  are  further  sorted  into  bands  of  1  m/sec, 
having  their  upper  and  lower  bounds  at  integral  wind  velocities. 
The  data  falling  into  each  band  are  first  averaged, and  tneir 
standard  deviations  from  the  average  value  is  subsequently  de¬ 
termined.  Along  with  the  average  value  of  each  band,  a  short 
vertical  line  is  drawn  to  indicate  the  standard  deviation  of  the 
data  from  the  average  value  for  each  band;  see  Figure  2.  The 
average  percentage  of  tne  standard  deviation  from  the  average 
value  is  only  22  percent.  A  discontinuity  in  the  wind-stress 
coefficient  is  shown  very  clearly  at  Uio  =  15  m/sec.  For  wind 
velocities  below  this  discontinuity,  the  wind-stress  coefficient 
is  approximated  very  well  by  the  proposed  correlation  curve,  and 
above  this  discontinuity  the  wind-stress  has  a  constant  value. 
The  mechanism  responsible  for  this  discontinuity  will  be  ex¬ 
plained  in  a  latter  section. 

FLOW  REGIMES  AND  SURFACE  ROUGHNESSES 
UNDER  DIFFERENT  SEA  CONDITIONS 

According  to  Schlichting  (1955)*  the  surface  roughness  can 
be  found  from  [ l] 

k  =  30  h  =  - — - r-  [  a] 

ext>  (h/C  2  ) 

y 

It  is  obvious  from  the  calculated  surface  r  u.grn*  ss-  r  otted  1:. 
Figure  3a,  that  there  are  two  different  regions,  .’.ami.  '.y,  t  he  r-  - 
gio.n  of  roughness  establishment  In  which  s,<r  fnc-  i  v..gn: 
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increases  with  wind  velocity,  and  the  region  of  established 
roughness  in  which  surface  roughness  reaches  a  saturated 
va lue . 

The  roughness  Reynolds  numbers  which  characterize  the 
boundary  layer  flow  conditions  can  now  be  obtained  from  the  cal¬ 
culated  roughnesses  and  are  plotted  in  Figure  3b.  In  the  same 
figure,  the  limits  for  different  flow  regimes,  (Schlichting, 
1955)  —  aerodynamica lly  smooth  flow,  ku*/v  <  5;  transition, 

5  <  ku*/v  <70;  aerodynamica lly  rough  flow,  ku*/v  >70  — 

are  also  marked,  where  v  is  the  kinematic  viscosity  of  air.  The 
boundary  layer  flow  conditions,  found  from  Figure  3b,  are  sum¬ 
marized  oelow. 


TABLE  1 

Boundary -Layer  Flow  Conditions 
at  Air-Sea  Interface 


Boundary  Layer  Flow  Condition 

Wind  Velocity 
Range 
m/sec 

Aerodynamica lly  Smooth  Flow 

Transition 

Aerodynamica lly  Rough  Flow 

U10  <  3 

3  <  U10  <  7 

U10  >  7 

Since  boundary  layer  conditions  and  surface  rougnnesses  are 
the  most  frequently  discussed  subjects  concerning  the  air-sea 
interface.  It  is  surprising  that  only  piecemeal  information 
containing  widely  scattered  values  have  been  reported.  The 


t 
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present  result  Is  believed  to  be  the  first  systematic  presen¬ 
tation  of  the  oceanic-determined  surface  roughness  and  subse¬ 
quently  the  boundary  layer  flow  conditions . 

The  roughness  at  the  air-sea  interface  is  composed  of  a 
great  variety  of  moving  waves,  differing  in  height,  shape  and 
phase  velocity.  Moreover,  these  waves  are  subjected  to  contin¬ 
uous  and  irregular  changes.  The  Instantaneous  picture  of  the 
complicated  sea-surface  structure  shows  short  waves  superim¬ 
posed  upon  longer  waves  of  various  sizes. 

If  we  consider  the  maximum  wave  height -t o-length  ratio  to 
be  0.14,  then  we  find  that  the  minimum  wave  lengths  correspond¬ 
ing  to  the  estimated  surface  roughnesses  in  the  aerody- 

namically  rough  regime  are  between  2.1  and  58  cm,  which  fall 
within  the  gravity  wave  category.  At  low  wind  velocities 
(3  m/sec  <  U10  <  7  m/sec),  the  minimum  wave  length  corresponding 
to  the  respective  surface  roughness  has  a  very  small  value.  How¬ 
ever,  it  is  clear  from  Schlichting  (1955)  that  the  height  of  thc- 
actual  roughness  is  usually  greater  (sometimes  much  greater)  than 
that  of  the  surface  roughness  calculated  from  the  wind  profile, 
depending  on  the  density  of  the  actual  roughness.  On  the  other 
hand,  an  Indirect  proof  of  relt  ting  the  surface  roughness  to 
gravity  waves  at  low  wind  velocities  is  provided  by  showing  the 
applicability  of  Cnarnbck's  nondimensiora 1  expression  to  the 
wind-stress  data  compiled  within  this  wi..d  velocity  range;  see 
Figure  2  and  discussions  in  lattr  sections.  The  restriction  of 
the  applicability  of  the  Charnock  expression  to  the  cases  of  grav¬ 
ity  waves  governing  surface  roughnesses  has  been,  r-stal- lisned  Lv 
Wu  (1967). 
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More  analytical  studies  along  with  detailed  oceanic  obser¬ 
vations  are  necessary  to  understand  the  precise  nature  of  the 
surface  roughness  and  its  variation  with  sea  conditions,  espe¬ 
cially  for  the  low  wind  velocity  range,  discussed  above,  and 
for  the  regime  of  established  roughness.  Tne  results  here  seem 
to  relate  the  surface  roughness  to  the  short  gravity  waves  and 
demonstrate  no  direct  correlation  between  the  surface  roughness 
and  the  mean  square  surface  slope;  see  Figure  4.  The  latter 
correlation  was  suggested  earlier  by  Munk  (1955)- 

REGIMES  OF  WIND-STRESS  AT  AIR-SEA  INTERFACE 

We  have  shown  the  boundary  layer  flow  conditions  at  the 
air-sea  interface  for  various  wind  velocities.  As  shown  in  ^Ig 
ure  2,  In  which  the  scales  of  sea  state  as  well  as  of  Beaufort 
number  are  marked,  the  proposed  correlation  formula  approxi¬ 
mates  the  trend  of  the  oceanic  data  at  the  frequently  occurring 
light  wind,  (3  m/sec  <  Ui  o  <  II  m/sec,  or  sea  state  1-4  corres¬ 
ponding  to  Beaufort  number  1-6).  For  a  wind  velocity  less  than 
3  m/sec  or  breeze,  (sea  state  0-1  and  Beaufort  number  0-1),  the 
flow  is  aerodynamics lly  smooth,  and  a  higher  wind-stress  coeffi 
clent  is  shown.  The  wind-stress  coefficient  remains  constant 
for  strong  wind,  (Uio  >  15  m/sec,  or  sea  state  above  5  corres¬ 
ponding  to  Beaufort  number  above  7).  A  discussion  of  each  case 
follows . 
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Wind -Stress  for  Breeze:  The  wind -stress  coefficient  for 
aerodynamically  smooth  flow  can  be  shown  as 


and  is  plotted  as  a  dotted  line  in  Figure  2.  As  Indicated  by 
the  compiled  data,  the  flow  is  aerodynamically  smooth  for 
Uj.o  3  m/sec  as  concluded  in  the  previous  section.  Tne  wind- 
stress  coefficient  in  the  transition  region  (3  m/sec  <  Ui o  <  7 
m/sec),  however,  is  snown  to  be  represented  very  well  by  the 
formula  for  aerodynamically  rough  flow  rather  than  (the  dotted 
line  corresponding  to)  aerodynamically  smooth  flow. 

Wind  Stress  for  Light  Wind:  In  Figure  2,  the  region  of 
applicability  of  the  Charnock  relationship  is  shown  clearly. 

The  proportionality  constant  of  the  Charnock  expression  was  de¬ 
termined  by  Wu  (1967)  in  a  laboratory  channel  as  0.0112  and  has 
been  shown  in  Figure  1  to  fit  the  field  wind-stress  data  fairly 
well.  However,  at  this  stage  it  may  be  worthwhile  to  determine 
a  new  constant  based  on  the  data  compiled  from  thirty  indepen¬ 
dent  oceanic  investigations. 

This  constant  is  determined  Indirectly  by  finding  a  vain- 
which,  after  substituting  Into  [3],  offers  tne  best  fit  'f  the 
stress  coefficient  curve  to  the  compiler)  data  ir.  th«  region  of 
light  wind.  Equal  weight  has  been  assigned  to  each  original 
data  point;  consequently,  c-ach  averaged  data  point  smwr.  in 
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Figure  2  possesses  a  different  weight  proportional  to  the  num¬ 
ber  of  the  original  data  points  it  contains.  The  root-mean- 
square  deviation  of  the  compiled  data  points  from  the  curves 
with  the  substitution  of  various  values  for  "a"  are  plotted  in 
Figure  5a.  The  minimum  point  corresponds  to  a  new  Charnock  pro¬ 
portionality  constant  determined  directly  from  oceanic  data 


a  =  O.OI56  =  — — 

us*/g 


[6] 


The  wind-stress  coefficient  for  this  most  frequently  encountered 
regime  can  now  be  written  as 


C 

y 


•bn 


_y£_ 


^0.0156  c  u  3 
y  y 


[7] 


and  is  drawn  as  a  continuous  curve  in  Figure  2;  in  which  an  ex¬ 
cellent  representation  of  the  compiled  data  by  this  expression 
is  shown. 


The  expression  for  the  wind  velocity  distribution  can  then 
be  obtained  as 


u#  * 


+  10.4 
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Wlnd  Stress  for  Strong  Wind :  Beyond  sea  state  5 
Ui  o  >  15  m/sec)  the  wind -stress  coefficient  shows  a  trend  to 
remain  constant,  which  is  obvious  from  [  l]  ,  suggesting  the 
saturation  of  roughness,  p.  More  studies,  especially  oceanic 
observations,  are  needed  in  order  to  thoroughly  understand 
this  mechanism.  However,  this  trend  does  not  conflict  with  ex¬ 
pectation  since  the  waves  responsible  for  the  surface  roughness 
cannot  grow  with  the  wind  velocity  forever.  As  the  wind  blows 
harder,  the  transfer  of  more  energy  from  the  air  merely  provides 
energy  to  waves  with  large  wave  lengths,  which  by  having  very 
high  phase  velocities  and  relatively  flat  shapes,  is  hardly  re¬ 
sponsible  for  the  surface  roughness. 

The  wind-stress  In  this  case  is  found  to  be 


C  =  2.6  x  10“ 3 

y 


l9l 


It  may  be  worthwhile  to  note  that  this  trend,  remaining  con¬ 
stant  at  high  wind  velocities,  is  also  shown  in  the  Nikuradse 
experiments  for  rough  pipe  flow  at  very  high  Reynolds  numbers. 
These  experiments  bear  a  close  analogy  with  the  present  case. 


DISCONTINUITY  OF  WIND -STRESS 
COEFFICIENT  —  CRITICAL  WIND  VELOCITY 


As  stated  in  previous  sections,  mere  oceanic  otn-Tva*.  i  nc 
of  the  microscopic  structure  of  the  air-sen  lnt  t  rl’ac  •  ni  •  •  >.<  •  i*.  d 
in  order  to  understand  the  precise  nature  of  t  h>  ,  .:T>e-  i  .t->  - 
ness  which  in  turn  governs  the  wind -stress .  VC*.  •*,  ,  .  •>*.  \  •  mi  t  i 
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to  employ  oceanic  Information  already  available  to  explain  the 
existence  of  the  discontinuity  of  the  wind-stress  coefficient 
at  Uio  =  15  m/sec. 

A  well  coordinated  and  frequently  quoted  presentation  of 
oceanic  data  concerning  the  wave  period,  phase  velocity,  and 
significant  wave  height  was  presented  by  Wlegel  (1962);  its 
small  scale  reproduction  is  shown  in  Figure  6.  Empirical  for¬ 
mulae  for  the  dimensionless  phase  velocity  (C/U10)  and  the  sig¬ 
nificant  wave  height  (gH^/U3 1 0 )  at  various  combinations  of  wind 

velocity  and  fetch  (gF/l^io)  are  found  to  be 


C 

Uio 


0,050 


gF 

U3  1  0 


0.30 


U*io 


0.0031 


U®io 


[  10] 


[11] 


Data  points  for  (gH^/U^io)  are  seen  to  deviate  from  the  drawn 
line  at  the  left  half  of  the  plot;  however,  the  oceanic  obser¬ 
vations  , owing  to  their  great  fetches,  occupy  only  the  right 
half  of  the  plot. 

From  these  two  empirical  formulae,  the  average  phase  ve¬ 
locity  and  the  significant  wave  height  for  various  wind  veloc¬ 
ities,  at  any  given  fetch,  can  be  obtained.  On  the  other  hand, 
based  on  the  logarithmic  velocity  distribution  and  the  compiled 
wind-stress  coefficient  data,  the  wind  velocity  at  a  height  of 
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the  significant  wave  amplitude,  (half  the  significant  wave 
height),  from  the  mean  water  level,  can  be  determined  from  C 8] . 

A  sample  plot  of  these  two  velocities,  the  average  phase  veloc¬ 
ity,  C,  and  the  wind  velocity  at  the  height  of  the  significant 
wave  amplitude  (Uio ).  for  a  fetch  of  50  km  is  shown  in  Fig- 

ure  6a.  Notice  that  these  two  velocities  intersect.  Let  us 
call  U  the  wind  velocity  at  which,  for  a  particular  fetch,  the 
intersection  occurs.  Following  tne  same  procedure,  the  inter¬ 
secting  wind  velocities  at  various  fetches  are  found  and 
plotted  in  Figure  6b.  Surprisingly  enough,  as  the  fetch  in¬ 
creases,  this  wind  velocity,  U^,  approaches  asymf t ot ica 1 ly  to 
Uio  =  15  m/sec.  This  is  the  very  velocity  where  the  discontin¬ 
uity  of  the  wind-stress  coefficient  occurs.  These  fetches 
(50  km  or  greater)  where  the  intersecting  velocity  approaches 
15  m/sec  are  very  common  for  oceanic  applications. 

The  implication  of  the  intersection  of  these  two  veloci¬ 
ties  (wind  velocity  and  phase  velocity)  is  that  the  typical 
waves  pull  the  air  mass  for  U10  less  than  15  m/sec  and  the  air 
mass  pushes  the  typical  waves  for  Ui  c  greater  than  15  n/sec.  A 
drastic  change  in  the  structure  of  the  air  flow  and,  cm.'oqui  r.t  l,v . 
in  the  wind-stress  coefficient  is  thus  inevitable.  This  finding 
provides  the  explanation  for  the  existence  of  the  disc. ntinol ty 
of  the  wind -stress  coefficient.  In  addition,  this  finding  ff  rs 
encouragement  tr  oceanographers  by  proving  the  cr.si.-t  i-rc.v  f  t  n- 
wind-stress  data  (this  discontinuity  was  sh  ■ -wn  a  Is;  •*.  in.  >  y.  ry 
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set  of  wind-stress  data)  and  the  consistency  between  the  wind- 
stress  data  and  other  oceanic  observations  (by  using  wave 
measurement  results  to  explain  the  trend  of  wind-stress  data). 

There  are  two  commonly  used,  but  very  confusing,  geophys¬ 
ical  terms:  light  wind  and  strong  wind.  Various  values  of 
wind  velocities,  between  8  and  15  m/sec,  are  quoted  as  the 
separating  velocity.  The  present  finding  thus  proposes  a  nat¬ 
ural  division  (Uj.0  =  15  m/sec)  for  the  light  and  strong  winds. 

This  finding  may  also  be  related  to  the  so-called  "critical 
wind  velocity"  first  proposed  by  Munk  (19^7)  but  accorded  er¬ 
roneously  to  the  transition  of  the  boundary  layer  flow  from 
laminar  to  turbulent. 

PROPOSED  WIND-STRESS  FORMULAE 

Many  reviewers  have  proposed  various  wind-stress  formulae. 

They  generally  contain  only  a  few  sets  of  data  which  may  not  be 
representative  and,  more  Importantly,  make  little  attempt  to  pro¬ 
vide  some  physical  reasoning  about  the  grouping  of  the  data.  Some 
Investigators  Improperly  compiled  the  laboratory  data  with  the  oce¬ 
anic  results  by  employing  y  =  10  cm  for  the  former  and  y  =  10  m 
for  the  latter.  Other  investigators  neglected  the  obvious  trend 
of  the  dependency  of  the  wind-stress  coefficient  on  wind  velocity. 
The  discontinuity  of  the  wind-stress  coefficient  was  not  noted  in 
most  of  the  reviews.  The  present  review  has  separated  the  oceanic 
observations  from  the  laboratory  results  (see  Figure  l)  and 
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preserved  the  dependency  of  tne  wind -stress  coefficient  on  the 
wind  velocity  (see  Figure  2).  Moreover,  the  finding  of  the 
boundary  layer  flow  conditions,  aerodynamica lly  smooth  or  rough 
(see  Figure  3),  and  of  the  discontinuity  of  the  wind-stress  co¬ 
efficient  (see  Figure  2)  help  very  much  in  the  proper  handling 
of  the  wind-stress  data.  Some  of  the  wind-stress  formulae  dis¬ 
cussed  in  the  previous  section  have  not  been  expressed  explic¬ 
itly  and  simply  for  common  applications.  Consequently,  the 
following  approximations  of  the  wind-stress  formula  are  sug¬ 
gested  for  various  cases  and  are  compared  with  some  available 
authoritative  formulae  at  the  end  of  this  section. 

Breeze :  Following  the  frictional  resistance  formula  for 

the  turbulent  boundary  layer  along  a  smooth  plate,  the  wind- 
stress  coefficient  for  this  case  can  be  approximated  by 

C10  =  x  10' 3  for  Uio  <  1  m/sec  Cl2] 

i/ 

Uio  4 

Light  Wind:  In  order  to  approximate  the  trend  of  the 
curve  expressed  in  [ 7 3  and  to  select  the  simplest  formula  pos¬ 
sible,  the  following  Is  suggested 

i 

Ci  o  =  0.5  x  Uioa  x  10' 3  for  1  m/sec  <  Uio  <  15  m/scc  C  13  3 

This  curve,  drawn  in  Figure  2,  Is  seen  to  approximate  very  well 
the  original  wind-stress  curve  except  at  the  lower  wind  veloc¬ 
ity  end.  As  shown  in  the  last  section,  the  boundary  layer  flow 
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condltion  Is  In  the  transition  region  at  this  lower  wind  veloc¬ 
ity  end.  The  uncertainty  of  the  data  in  the  transition  region 
makes  this  slight  lation  even  more  discernable. 

Strong  Wind:  The  wind-stress  coefficient  is  already  in  a 
very  simple  form  as  shown  in  [93.  This  formula  is  applicable 
for  a  wind  velocity  greater  than  15  m/src. 

The  present  and  the  previously  proposed  and  often  quoted 
formulae  are  tabulated  in  the  following  table.  The  root-mean- 
square  errors  of  the  wind-stress  data  from  the  proposed  curves 
are  also  presented  for  the  purpose  of  comparison.  In  the  error 
calculation,  equal  weight  once  again  has  been  assigned  to  every 
original  data  point. 


TABLE  2 

Comparison  of  Wind -Stress  Formulae 


Author 

Year 

Proposed  Formula 

Ci  0  x  103  = 

Range  of 
Application 
m/sec 

Root -Mean- 
Square  Error 

Wu 

1967 

1 

0.5  x  U10  2 

1  <  Uxo  <  15 

0. 103 

2.6 

U10  >  15 

0. 130 

Deacon  & 
Webb 

1962 

1.0  ■+•  0 . 07  x  Ui  0 

-Tf- 

1 — 1 

v 

0 

tS 

V 

0.246 

Wi Ison 

I960 

1.49 

1  <  U10  <  10 

0.369 

2.37 

U10  >  10 

0.403 

Sheppard 

1958 

0.8  +  0.114  x  U10 

0 . 3  67 

Francis 

1951 

1*3  x  Ui  0 

bdessei 

0.379 

HHMB 

1948 

1 

0.9  X  Ul0 

1  <  U10  <  :o 

1.441 
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CONCLUSION 

The  wind-stress  data  are  compiled  from  the  following  ind 
pendent  laboratory  and  field  investigations 


TABLE  3 

Sources  of  Wind  Stress  Data 


Laboratory  Investigations 

rield  Investiga 

t  ions 

Author 

Year 

Author 

Y  ear 

Francis 

1951 

Colding 

1876 

Francis 

1951 

Ekman 

19C5 

Keulegan 

1951 

J  ohnson 

1927 

Johnson,  et  al 

1952 

Shoulejkin 

1918 

Hamada 

1953 

Montgomery 

1936 

Sibul 

1953 

Sutcliffe 

1936 

Hayami,  et  al 

1959 

Palmen,  et  al 

1938 

Fitzgerald 

1963 

Bruch 

1990 

Fitzgerald 

1963 

Scha lkwijk 

1997 

Fltzgera Id 

1963 

Hela 

1998 

Hidy,  et  al 

1966 

Roll 

1998 

Wu 

1967 

Corkan 

193v 

Durst 

1959 

Sheppard 

1 951 

Sheppard,  et  al 

1952 

van  Dorn 

1952 

Hel 1st rom 

1°33 

Hell strom 

1953 

Keulegan 

i  -35 

Hay 

1  -.55 , 

Mcllroy 

1 

Cnarnock,  et  al 

1  'Re 

Deacon,  et  a  ! 

I  ‘Re 

Ta  k  a  ha  sh  i 

1  158 

Frocks 

1  >5° 

Vinogradova 

1  15a 

Bruce,  et  a! 

1  h.  ; 

Brocks 

]  9b. 

Deacon 

1  H 

Deacon 

1  V 

he  results  found  in  this  report 


a  re 


mma riued 


HYDRONAUTICS,  INCORPORATED 


-20- 


■*r 

LU 


U 

< 


Z 

< 

LU 

UO 

I 


APPROXIMATE  WIND  - 
STRESS  FORMULA 
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FIGURE  1  -  SUMMARY  OF  WIND-STRESS  MEASUREMENTS  IN  LABORATORIES  AND  IN  FIELD 
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FIGURE  2  -  OCEANIC  MEASUREMENTS  OF  WIND-STRESS  UNDER  DIFFERENT  SEA  CONDITIONS 


ROUGHNESS  REYNOLDS  NUMBER,  —  SURFACE  ROUGHNESS , 


HYDRONAUTICS,  INCORPORATED 


E 

u 


* 


WIND  VELOCITY,  Ujq,  m/sec 


;a)  SURFACE  ROUGHNESS  AT  AIR-SEA  INTERFACE 


WIND  VELOCITY,  U,0,  m/sec 


(b)  BOUNDARY  LAYER  FLOW  CONDITIONS 


FIGURE  3  -  SURFACE  ROUGHNESS  AND  FLOW  REGIMES  UNDER  DIFFERENT 
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FIGURE  4  -  SURFACE  ROUGHNESS  AND  MEAN  SQUARE  SURFACE  SLOPE  UNDER 
DIFFERENT  CONDITIONS  ( THE  MEAN  SQUARE  SURFACE  SLOPES 
ARE  REPRODUCED  FROM  COX  AND  MUNK,  1956) 
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FIGURE  6  -  CHANGING  OF  WIND-STRESS  REGIMES 
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